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Abstract 

The one-dimensional (ID) t — J model is investigated using the density ma- 
trix renormalization group (DMRG) method. We report for the first time a 
generalization of the DMRG method to the case of arbitrary band filling and 
prove a theorem with respect to the reduced density matrix that accelerates 
the numerical computation. Lastly, using the extended DMRG method, we 
present the ground state electron momentum distribution, spin and charge 
correlation functions. The 'ikp anomaly of the momentum distribution func- 
tion first discussed by Ogata and Shiba is shown to disappear as J increases. 
We also argue that there exists a density-independent Jc beyond which the 

system becomes an electron solid. 
71.10.-hx,71.27.-ha,75.10.Jm 
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I. INTRODUCTION 

For a number of years, the more so since the discovery of high-T^ superconductors!, the 
study of strongly correlated electron systems has been a major theme of theoretical con- 
densed matter physics. The t — J model is one of the simplest of such model^. Although 
high-Tc cuprates are at least two-dimensional {2D) systems, it is also relevant to fully un- 
derstand the one-dimensional (1-D) model. The Hamiltonian for the t — J model in one 
dimension can be written in the subspace of no doubly occupied sites as 

H = -tJ2{cLci+ia + cf+i^Q^) + J ^(SrSi+i - ^riiUi+i), (1) 

ia i 

where and Qo- are the creation and annihilation operators for an electron at lattice site i 
with spin a, while Sj and nj are the corresponding spin and electron number operators. This 
model has been solved exactly only for J — 0, where it is equivalent to the U oo Hubbard 
model, and at the supersymmetric point! J = 2t. In both cases the ground state at arbitrary 
density belongs to a broad class of interacting Fermi systems known as Luttinger liquids!, 
which exhibit power-law decay of correlation functions with exponents characterized by a sin- 
gle parameter!. Additionally, for very large J/t, the attractive Heisenberg interaction term 
dominates the kinetic energy and the model phase separatesl'i. The level of understanding 
of this model is derived mostly from small cluster exact diagonalizationsi, variational Monte 
Carlo methods! and finite temperature Monte Carlo simulations for relatively larger system 
sizes!. There still lurks, however, the question of whether the thermodynamic limit has been 
reached or not. In this paper we use the density matrix renormalization group (DMRG) 
method0 to study the ground-state properties of the t — J model in the thermodynamic 
limit. 

The text is organized as follows. In Sec. II, we review the formulation of DMRG 
and its extension to arbitrary band-filling . In Sec. Ill, we discuss our numerical results. 
We present ground-state static spin and charge correlation functions as well as electron 
momentum distribution functions. Finally in Sec. IV we summarize our results. 
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II. DMRG FORMULATIONS 

The DMRG technique was developed by White0 in 1992. It leads to highly accurate 
results for systems much larger than those which can be solved by exact diagonalization. The 
DMRG allows a systematic reduction of the Hilbert space to the basis states most relevant to 
describe a given eigenstate {e.g. the ground state) of a large system. This is to be contrasted 
with previous real space renormalization techniques in which the lowest states are kept. A 
general iteration step of the method for open boundary conditions proceeds as follows: a) 
The effective Hamiltonian defined for the superblock 1 + 2 + 2' + 1' (where the block 1 and 1' 
come from previous iterations and 2 and 2' are new added ones) is diagonalized to obtain the 
ground-state wave function \ip > (other states could be also kept), b) The reduced density 
matrix of blocks 1 + 2: pj j/ = J2j i^iji^i'j is constructed, where ipij =< i ® >, the states 
\i > {\j >) belongs to the Hilbert space of blocks 1 and 2 (1' and 2'). The eigenstates of 
p with the highest eigenvalues (equivalent to the most probable states of blocks 1 + 2 in 
the ground state of the superblock) are kept up to a certain cutoff, c) These states form a 
new reduced basis in which all the operators have to be expanded and the block 1 + 2 is 
renamed as block 1. d) A new block 2 is added (one site in our case) and the new superblock 
(1 + 2 + 2' + 1') is formed as the direct product of the states of all the blocks (the blocks 
1' and 2' are identical to blocks 1 and 2 respectively). The method has been successfully 
applied to problems such as the Haldane gap of spin-1 chains, critical exponents of spin-i 
chains, the ID Kondo-insulator and two-chain Hubbard mo del0. 

In our study of the ID t — J model we have used the infinite-size version of the above 
iteration scheme to reach the thermodynamic limit. One immediately realizes, however, there 
is a problem in keeping the electron density fixed in the iteration process since we insert 
only two sites at each iteration (which only makes the half-filling and quarter-filling cases 
invariant). To get around this, we construct the reduced density matrix from two ground 
states that bracket the desired density. To be more specific, if our desired electron density is 
n and the current superblock lattice size is iV, then we can always find two nearest integers 
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A^i and N2 such that A^i < nN < N2. Assuming \'il){Ni) > is the ground state wavefunction 
of Ni electrons and |'?/'(A^2) > is that of A^2 electrons, we build the reduced density matrix 
by the following weighting procedures: 

Ay = Pi ^.,(iVi)^;,^.(iVi) + p2 ^.,(iV2)^*^.(iV2) 
< nAT = pliVi + P2iV2 (2) 
1 = Pi + P2 

It is clear that the above construction ensures the desired constant band filling at every 
iteration. The success of this construction inherently requires that the ground state be 
homogeneous. Therefore one should expect failure when the system goes into the phase 
separation regime. However, phase separation can still be studied using the finite size 
version of DMRG0 which does not require translational invariance. It is perhaps worthwhile 
to remark that in our computer program we have generally targeted three states since the 
ground state has two-fold degeneracy when the electron number is odd. 

The correlation functions are calculated when the iteration has converged. Typically we 
start to measure them after the superblock size reaches 40 lattice sites. More specifically we 
measure the correlation between operators in the middle of the superblock. For example, 
the nearest-neighbor correlation is measured between sites 20 and 21 in the superblock size 
40; while the second-nearest-neighbor correlation is measured between sites 20 and 22 in the 
next iteration (i.e when superblock size becomes 42). Similarly the third- nearest-neighbor 
correlation is measured between sites 20 and 23 in the superblock size 44. Thus as we iterate 
our procedure to increase the system size we obtain correlations over increasing distances. 
Finally, the desired ground-state correlations are obtained in a similar way as the reduced 
density matrix, e.g. the spin-spin correlation is measured in the following way 

^(r - r') = pi < V^(A^i)|S(r) ■ > +p2 < V^(Ar2)|S(r) ■ Sir'^N^) > . (3) 

Before we discuss the numerical results, let us prove the following theorem: the eigen- 
states of the reduced density matrix retain their good quantum numbers when the correspond- 
ing operator of the superblock is a direct sum of the operators of its subsystems. The electron 
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number operator N and the total z-component of the spin operator Sz are two such oper- 
ators. To prove such a theorem all we need to show is that the reduced density matrix is 
block diagonal in its good quantum operator subspace. Let us consider one subspace at a 
time, e.g. the electron number subspace. Then the most general form of the eigenstate of 
the superblock with L electrons can be written as 

|V^(L)>= Ec(^i'«;^2,/3)|/i,a>i®|/2,/?>2, (4) 

h+l2=L a,p 

where >i stands for the basis wavefunction of block 1 with li electrons, and a labels 
other quantum numbers, while \l2il3 >2 represents the basis wavefunction of block 2 with I2 
electrons, and (3 marks other quantum numbers. Thus the reduced density matrix elements 
have the form 

= c(/i, «; P) " c*{l[, a'] h, P). (5) 

12,I3 

However according to the construction of the wavefunction we have both L = li + I2 and 
L = l[ + I2 leading automatically to the conclusion li = l[ . Therefore the reduced density 
matrix is block diagonal in the electron number subspace. In the same way one can show 
that the reduced density matrix is block diagonal in the Sz subspace. We emphasize that 
this theorem can be explicitly implemented to accelerate the diagonalization of the reduced 
density matrix in the subspace of its good quantum numbers. 

III. NUMERICAL RESULTS 

For most of the numerical results reported here we have kept 110 stateslil in blocks 1 and 
1', the truncation error defined as 1 —p{m) (where p{m) is the summation of the highest m 
eigenvalues of the reduced density matrix) is of the order of 10~^. In Table I we list some 
of the ground-state energies per site Eg as a function of band filling n and spin exchange 
coupling J. In the same table we also list the results obtained from the exact Bethe ansatz 
solution at J = 2t. As can be seen immediately from the table, our DMRG results are 
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highly accurate. To make it convenient for the reader we quote the coupled Bethe ansatz 
equationsll determining the ground-state energy at density n and J = 2t: 

Q 

J p{u)di' = 1 — n (6) 
-Q 

Q 

p{u) = 2R{2u) + j 2R{2[u - u'])p{u')dv\ (7) 

-Q 

where R{x) denotes Shiba's function 

R[x) = — / doj- ^. (8) 

— oo 

Then the ground-state energy is given by 

Eg = 2t[l-n-7rp(0)]. (9) 

Since we are going to present our results in momentum space it is now opportune to 
discuss the way we analyse our data. First, because open boundary conditions are used 
in our DMRG procedure, we obviously loose translational invariance. To overcome this 
boundary effect we have done an average on the real space correlation functions before we do 
the Fourier transform into momentum space. More specifically we obtain several real-space 
correlation functions of the same system by starting the measurement at different superblock 
sizes. For example we start our real-space density- density correlation pp{r) =< riiTii+r > 
measurement after the superblocks reach the sizes of 40, 42, 44, 46 and 48 sites and then 
take the average among them. The final Fourier transform is done in the usual way: 

-I N N 

PP^k) = ^ E E e*'('-"H< >-<ni><n, >], (10) 
1=1 h=i 

where N is the largest separation available in our measurement. For the results of this work 
we stop at a superblock size of 150 sites. 

In Fig. 1(a) we plot the electron momentum distribution function for quarter filled band 
at J = O.lt. We see that this result is almost identical to the ID infinite-t/ Hubbard model 
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studied by Ogata and Shiba£3. This is not surprising considering the fact that the small J 
limit of the t — J model is equivalent to the strong coupling limit of the Hubbard model. Here 
we have both kp (= wn j'i) and Zkp anormalies as indicated by arrows in the figure. Unlike 
normal Fermi liquids, however, the ?)kp anomaly is a new feature which is related to the effect 
of the coupled holon and spinon representation of the normal electronlll. Fig. 1(b) shows the 
spin-spin correlation. There is a anomaly refiecting the nature of the antiferromagnetic 
exchange among nearest-neighbor electrons. Fig. 1(c) is the density- density correlation 
function. Here one observes the ^kp (= vr) anomaly only. 

Before discussing the effect of increasing J let us remember the two competing factors 
in correlation effects: a) The nearest-neighbor antiferromagnetic spin exchange J favors 
the formation of spin-singlet electron pairs on nearest-neighbor sites which in turn leads to 
2kp spin density wave (SDW) fluctuations. But eventually, when J is large enough, there 
emerges a phase separation exhibiting antiferromagnetic SDW fluctuations. 6) The kinetic 
energy t term favors the delocalization of the electrons. But the Pauli principle as well 
as on-site repulsion promote uniform separation between electrons synonymous with Akp 
charge density wave (CDW) fluctuations. On the other hand the formation of bound singlet 
pairs due to J term enhances Ikp CDW fluctuations. 

In Fig. 2(a) it is seen that as the antiferromagnetic exchange J increases, the momentum 
distribution function is modifled more drastically in the region k > kp where it changes from 
decreasing to increasing as a function of k. Notably the 3kp feature has been washed out 
for larger J. Besides, the critical exponent! at the Fermi wavevector kp seems to decrease 
as J increases. In Fig. 2(b) we have plotted the spin-spin correlation function at quarter 
fllling with J/t = 0.5,2, and 2.5. We see that at J = 0.5t there is a clear 2kF anomaly, 
however at J = 2t this anomaly is signiflcantly weakened. At J = 2.5t the 2kp anomaly 
seems to have completely disappeared and the maximum is located at = vr = k^p. This 
is an indication of singlet pairing between nearest-neighbor electrons. Fig. 2(c) shows 
the corresponding density- density correlation. At J = 0.5t there is only a Akp anomaly 
which indicates the system is quite uniformly distributed. However as one increases J the 
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Akp feature disappears and a new feature at 2kp develops which reflects the formation of 
nearest-neighbor pair. At J = 2.5t the behavior near k = has a tendency to flatten out. 
But the 2kp feature still dominates which shows singlet pairing. Fig. 3 shows similar results 
for the band filling n = 0.8. Again it is seen from Fig. 3(a) that J has the largest effect on 
the momentum distribution function in the regime k > kp. Fig. 3(b) is very much the same 
as of Fig. 2(b): as J increases the maximum in the spin-spin correlation function moves 
to the antiferromagnetic wavevector kAp = tt. But in Fig. 3(c) the J = 3t density-density 
correlation shows a clearer tendency towards phase separation near A; = 0. Both 2k p and 
Akp anomalies are seen in the density correlation at intermediate J as predicted!. 

As mentioned in the previous section our DMRG procedure is not applicable to the 
study of the phase separated regime. But we know from other numerical calculations there 
exists a density — dependent J'^ indicating the onset of a phase separation. Here we wish to 
argue that there exists another critical Jc (that is larger than the max[J'^]) beyond which the 
system forms an electron solid, i.e., a phase with no hole inside but having a thin interface at 
which the boundary electrons can still evaporate into the vacuum (hole). This Jc should be 
independent of the band filling since once the electron solid is reached there is only a single 
phase boundary between the hole region and the solid, and that is density independent. By 
definition, the first instability of such an electron solid comes from the ability to dissolve a 
single hole inside. Accordingly can be estimated from finite size exact diagonalizations 
by comparing the ground-state energy of a 2A^-site Heisenberg chain with that of 2 holes 
in the 2iV-site t — J model. This has been done by Ogata et alB on a 16-site system and 
by Assaad et all! with path-integral Monte Carlo simulation. Both results seem to point 
towards Jc ~ 3.6t. Further, Yokoyama et all have hinted the existence of such a Jc in their 
recent variational Monte Carlo study oit — J model. The same thing is also implied in the 
work of Ammon el al0. 
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IV. SUMMARY 



We have extended the DMRG method to allow calculations at arbitrary band fiUing and 
proven a theorem regarding the reduced density matrix that could accelerate the numerical 
calculation. We then used the extended DMRG procedure to study the ID t — J model. 
It is found that our procedure gives highly accurate ground-state energy and correlation 
functions. It is shown that the "ikp anomaly first discussed by Ogata and Shiba in the 
electron momentum distribution disappears as J is increased. It is argued that there exists 
a density-independent Jc beyond which the system forms an electron solid. 
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TABLES 

TABLE I. Ground-state energies of the IDt — J model calculated using DMRG method. Exact 

results for the J = 2t case are listed as a comparison. The last number in the brackets are the 
estimated value accurate at that digit. 



< n> 


J/t 


DMRG Eg/t 


Exact Eg/t 


0.4 


0.5 


-0.631(8) 




0.4 


1.0 


-0.664(2) 




0.4 


2.5 


-0.803(3) 




0.5 


0.1 


-0.647(0) 




0.5 


0.5 


-0.692(0) 




0.5 


2.0 


-0.903(6) 


-0.9036(4) 


0.5 


2.5 


-0.988(5) 




0.8 


0.1 


-0.116(2) 




0.8 


0.5 


-0.586(7) 




0.8 


2.0 


-1.246(4) 


-1.2464(4) 


0.8 


3.0 


-1.698(7) 




0.9 


1.0 


-0.756(4) 




0.9 


2.0 


-1.322(3) 


-1.3223(7) 


0.9 


3.0 


-1.891(9) 





10 



FIGURES 

FIG. 1. The ground state electron momentum distribution function (a), spin-spin correlation 
function (b), and density-density correlation (c) are plotted for quarter band filling at J = O.lt. 

FIG. 2. Same as that of Fig. 1 but with J = 0.5t, 2t and 2.5i. 

FIG. 3. Same as that of Fig. 2 with J = OM, 2t and 3t but at band filling 0.8. 
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